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Abstract. We analize the efficiency of wavefront reconstruction in the MultiConjugate 
Adaptive Optics system for the European Solar Telescope (EST). We present preliminary 
results derived from numerical simulations. We study a 4 meter class telescope with multiple 
deformable mirrors conjugated at variable heights. Along with common issues, difficulties 
peculiar to the solar case have to be considered, such as the low contrast and extended nature 
of the natural guide features. Our findings identify basic requirements for the EST Adaptive 
Optics system and show some of its capabilities. 
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1. Introduction 

The European Solar Telescope (EST), cur- 
rently in its design study phase, will be a 4- 
meter class solar telescope. EST will be opti- 
mized for the study of the magnetic coupling 
between the deep photosphere and the upper 
chromosphere. The main observational focus 
will be multi-wavelenght (spectro-)polarimetry 
at photospheric and chromospheric layers. 
Plasma phenomena on the solar scene seem to 
happen on a small scale (< 100km) and to ex- 
hibit high temporal evolution (seconds, min- 
utes), thus high spatial and temporal resolu- 
tions are needed. 

To fit these needs, EST will be equipped with 
a high technology Multi Conjugate Adaptive 
Optics (MCAO) system able to deliver high 
resolution images of the Sun with a cor- 
rected Field of View (FoV) of 1 arcm in 
dEST-Science Requirements docum ent 2008). 
It this paper we will report the first results of 
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our MCA O simulations, using a layer oriented 
approach (IRagazzoni et al. 20021) . 
In particular, we have studied how the cor- 
rected FoV depends on the configuration of the 
system tracking points on the solar surface and 
furthermore how it depends on the modal base 
dimension employed in the wavefront recon- 
struction. 

Finally we have tested through simulations a 
time delay reduction method based on wave- 
front forecasting. Its use seems to increase the 
performance of the MCAO system and thus the 
quality of the correction. 

2. MCAO simulations 

In this work we have analized the MCAO sys- 
tem behaviour using the LOST simulation tool 
dArcidiacono et al. 20041) along with some op- 
timization to include needs typical of the solar 
case and wavefront prediction capabilities. 
As mentioned above, we have used simulations 
to study the corrected FoV and the quality of 
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Fig. 1. SR maps for three different asterism radius: 7 (left), 10 (center) and 15 (right) arcsec. The 
best performance is achieved for the 10 arcsec radius for which the Strehl Ratio reaches 0.5 over 
a fairly large FoV (20 arcsec) and a value of 0.3 over 50 arcsec. In the right panel it is possible 
to see how further increasing the asterism radius produces an unsatisfactory correction in terms 
of FoV coverage and SR peak. Map FoV is equal to 1 arcmin x 1 arcmin 



the correction itself. 

To do this, we have simulated a MCAO system 
with the following characteristics: 

- 2 Deformable Mirrors (DM) 

- Conjugation altitudes: Km, 10 Km 

- Time step: 0.5 ms 

- Loop delay: 1 ms 

- Working wavelength: 555 nm 

- 7 Shack-Hartmann (SH) wavefront sensors 
(WFS) with 30x30 subpupils sampling 

- 2 Kolmogorov turbulence layers (FFT): 
km, 10 km 

- D/r Q : 26, 8 

where D is the telescope diameter and ro the 
Fried parameter. 

In these simulations we have used a hexagonal 
tracking point geometry (asterism). 
We have studied, at first, sky coverage depen- 
dence on asterism radius, using as a quality de- 
scriptor the Strehl Ratio parameter (SR): 



RS R dif fraction limited peak 

In fig. 1 we present SR maps for three dif- 
ferent choices of the asterism radius: 7, 10 and 
15 arcsec. 

As is clear, the best correction is achieved 
when using a 10 arcsec asterism radius, where 



the largest sky coverage is reached with a 20 
arcsec FoV and SR = 0.3. 
A SR — 0.3 is generally considered, in terms 
of spatial resolution, as a fairly good condi- 
tion where the PSF width is very close to the 
diffraction limited one. 

Every attempt to enlarge the corrected FoV by 
increasing the asterism radius produces an un- 
satisfactory condition with a SR that is not uni- 
form over the field and with a decrease of the 
SR peak itself (see right panel of fig. 1). Thus 
there exists a maximum useful asterism radius, 
above which the correction becomes unstable. 
Fig. 2 and 3 show, in fact, the maximum SR 
reached in every simulation iteration (time step 
= 0.5 ms) over the whole FoV for a 10 and 15 
arcsec asterism radius respectively. While in 
the first case the correction is stable, in the sec- 
ond case it is unstable: after 200 loops the SR 
peak decreases rapidly. 

The second case also shows a lower SR 
mean value with respect to the case in which 
the asterism radius is 10 arcsec. 

3. Time delay error reduction 

One of the most important problems occurring 
in the MCAO correction is the time delay 
error. 

Every MCAO system takes some milliseconds 
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Fig. 2. SR peak as a function of time for a 
10 arcsec asterism radius. Every iteration is 
equal to 0.5 ms. The correction is stable show- 
ing a fairly constant SR peak value. 
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Fig. 3. SR peak as a function of time for a 
15 arcsec asterism radius. Every iteration is 
equal to 0.5 ms. The correction is poorly sta- 
ble and after 200 loops the SR dramatically de- 
creases. 



to analize the incoming wavefront and com- 
pute the command matrix to move the mirror. 
This means that the correction doesn't exactly 
match the incoming wavefront, leaving some 
residuals which can decrease the SR and the 
sky coverage. 

One attempt to solve this problem is to use 
short-term wavefront forecasting. We have 
implemented a prediction tool algorithm based 
on linear extrapolation inside LOST. This 
algorithm is able to achieve a prediction of 
the modal coefficients starting from their 
respective known time series. 
In this way we have been able to simulate the 
behaviour of the prediction tool during the 




Fig. 4. SR plot as a function of time step. 
Comparison of SR without (red asterisks) and 
with (black boxes) linear prediction. When the 
linear predictor is operative the SR reaches a 
value of about 0.4 which is double the value of 
the case without prediction. 

loop closure. 

The algorithm makes use of the last two modal 
coeffiecient values to estimate the short-term 
predicted one. 

We have investigated the behaviour of the SR 
with and without the prediction tool with the 
same simulations parameters. The results are 
shown in fig. 4. 

As is clear, the SR is enhanced by the pre- 
diction, reaching a mean value of 0.4. In the 
same plot we show, for comparison, the perfor- 
mance of the correction without the linear pre- 
diction of modal coefficients: in this case the 
SR reaches 0.2 when the loop is closed. 

4. Modal base optimization using 
information theory 

Another way to minimize the time delay effects 
is by otpimizing the modal base. If it is possi- 
ble to use fewer modes in the correction then 
the calculation will be speeded up, leading to a 
reduction in the time delay. 
In this work we have optimized the standard 
Zernike modal base using information theory, 
working on real data obtained at the Ge rman 
VTT telescope ( von der Liihe et al. 20031) . 
We have reorder ed the moda l base by mu- 
tual information (lHavkin 20001) . forcing those 
modes which bring more information about the 
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Fig. 5. Fitting error as a function of the num- 
ber of modes used in the reconstruction. 
Mutual information Zernike reordering (con- 
nected points) leads to a faster decreasing trend 
with respect to standard modal ordering (aster- 
isk). For a given fitting error we can use far 
fewer modes in the case of mutual information 
ordering. 

incoming wavefront to appear in the first po- 
sitions. We have then estimated the fitting er- 
ror comparing it to that obtained in the case of 
standard ordering. 

The results of this analysis are shown in fig. 
5. Mutual information ordering leads to a faster 
decrease of the fitting error when compared to 
the standard Zernike ordering. In this way, for 
a given fitting error, we can choose to use fewer 
modes, thus speeding up the calculation. 
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